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ABSTRACT 

The proximity effect (PE) of QSOs is beheved to be useful in constraining the QSO lifetime. Obser- 
vations on the PE so far, however, give apparently contradictory results - some are consistent with a 
long QSO lifetime (> a few 10^ yr), but others appear to be only consistent with a short QSO lifetime 
< 10^ yr. In this paper, we show that this apparent contradiction may be solved by simultaneously 
taking into account both the effect due to the density enhancement in the QSO near zones and that 
due to the obscuration of the tori associated with the QSOs, using a large number of Monte-Carlo 
generated synthetic Lya forest spectra. We demonstrate that the QSO properties and environment 
can be constrained simultaneously by the transverse PE and the line of sight PE of bright type 1 QSOs 
together. The current available measurements on the PEs of type 1 QSOs suggest that (1) the density 
is significantly enhanced in the vicinity of the QSOs; (2) the QSO lifetime is consistent with being as 
large as a few 10^ yr and a substantially shorter lifetime (e.g., < 10^ yr) is not required; and (3) the 
half opening angle of the tori associated with QSOs is ~ 60°, consistent with some other independent 
estimates. Our simulations also show that the TPE of type 2 QSOs can be significantly different from 
that of type 1 QSOs, which may be useful to put further constraints on the QSO properties and the 
QSO environment. 

Subject headings: intergalactic medium-quasars: absorption lines-quasars: general-black hole physics 



1. INTRODUCTION 

The enormous UV radiation from a QSO can signif- 
icantly alter the ionization state of its surrounding in- 
tergalactic medium (IGM). The Lya absorption clouds 
near the QSO are more ionized than average, and the 
optical depth of Lya photons due to the absorbers 
was expected to decrease as the absorption redshift ap- 
proaches that of the QSO, which is named as the "prox- 
imity effect" (PE) (iCarswelll [19871: iBaitlik et aLlllQSl . 
The significance of the PE of a QSO is mainly de- 
termined by the following factors: (1) the strength of 
the metagalactic UV ionizing background (UVB); (2) 
the intrinsic properties of the QSO, including its lu- 
minosity, its age (or its detailed luminosity evolution), 
and the anisotropic feature in its UV radiation; and 
(3) the overdense environment of the IGM surround- 
ing the QSO. Therefore, the PE is useful n ot only to 
measure the UVB (e.g., iBajtlik et aH Il988l ). but also 
to constrain the lifetime/luminosity evolution of QSOs 
and the d ensity enhanceme nt in its near zones (e.g. , 

Jakobsenet al. 2003; Crofti l2003 iSchirber et all I200I 
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Adelberger 2004; Rollinde et al. 2005: Guimarae s et al.l 
2007; Kirkman & Tvtler 2008; Goncalves et al. 200^! 

The PE along the line of sight to a QSO (LOSPE) 
has been clearly detected through the high-resolution 
Lya absorption spectra of a combined QSO sample 
or some individual QSOs, which is wide ly used to 
estiruate the strength of t he UVB (e.g., IScott et all 
120001 iDairAglio et al.l l2008| ) and also the density en- 
hancement of the IGM surrounding the QSO (e.g., 
iRoUinde et al.ll2005HGuimaraes et al.ll2007l). For exam- 
ple. iDall'Aglio et al.l (|2008f) and iFaucher-Giguere et al.l 
(|2008b|) find that the UVB measured from the LOSPE 
(after including the effect of the overdense environment of 



QSOs) is consistent with those estimates obtained from 
other methods. The LOSPE is affected directly by the 
observed QSO luminosity, but almost irrelevant to the 
QSO's past luminosity evolution (including luminosity 
variation) and the age of its nuclear activity. In principle, 
the PE may be also affected by fluctuations in the QSO 
luminosity on short timescales ^ 2 x lO^yr, over which 
the perturbed ionization state of the highly ionized IGM 
can get back to the photoionization equilibrium. How- 
ever, the QSO variability studies give the average lifetime 
of an individual episode of QSO activity > lO^yr (e.g., 
iMartini fc Schneidei][2(ffl3t lMartinil[200l . 

The PE in the region other than along the line of 
sight to the QSO (i.e., the transverse PE or TPE) can 
be measured through the spectra of background QSOs 
(bgQSOs) or other bright UV sources whose light passed 
through the highly ionized region surrounding the fore- 
ground QSOs (fgQSOs). Unlike the LOSPE, the TPE 
should be sensitive to the age of the fgQSO, as a posi- 
tive detection of this effect requires that the light travel 
time from the fgQSO to the line of sight (LOS) to the 
bgQSO is smaller than the age of the fgQSO. Hence, the 
TPE has long been thou ght to be help ful in constr aining 
the QSO lifetime (e.g., Baitlik et al.' TgSS'; 'c7of?'200l 
Schirb er et al. 2004; Adelberger 2004; Kirkman & Tytlei 
I2OO8D . However, so far the results from observations on 
the significance of the TPE appear ambiguous. For ex- 
ample, many authors tried to search for the reduction 
in the optical depth or the number of HI Lya forest 
lines caused by the TPE through QSO pairs with small 
transverse separation (typically < a few Mpc), but they 
found no evidence of this reducti on a nd sug gested that 
the QSO lifetime < lO^yr (e.g., iFer nandez -Soto et al 
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120041: iKirkman fc Tvtleil I2008D . However, several au- 
thors have reported the detections of the TPE through 
the observations on the He II, rather than H I, Lya 
forest or metal absorption lines (e.g., Fjakobsc n et alj 
2Q03; 'Worscck & Wisotzki] 120061 : IWorseck et all 120071 : 
Goncalvcs et al. 2QQ8) . Those detections suggest that 



the QSO lifetime should be > lO^yr. Currently it is still 
unclear whether that apparent contradiction is due to 
some selection effects in observations or the complica- 
tions introduced by other factors, such as, the density 
enhancement in the QSO near zones or the anisotropic 
radiation from QSOs. 

In this paper, we study the LOSPE and the TPE of 
QSOs, with taking into account the density enhancement 
in the QSO proximity region, their anisotropic UV radi- 
ation, and their lifetime; and the purpose is not only to 
solve the apparent contradiction among the observational 
constraints on the QSO lifetime from current TPE ob- 
servations, but also to demonstrate that the LOSPE and 
the TPE combine to put constraints on both the QSO 
lifetime and the opening angle of torus associated with 
QSOs simultaneously. The paper is organized as follows. 
In Section [51 we illustrate the geometry of the proxim- 
ity region of a QSO in the observer's rest frame that 
can be affected by UV photons emitted from the QSO. 
Accounting the PE of (fg)QSOs, we use a Monte-Carlo 
method to generate a large number of Lya forest spectra 
of fgQSOs and bgQSOs as described in Section [31 based 
on some statistical distributions of the HI column den- 
sity, the Doppler factor of Lya absorption lines, and the 
reds hift distribution of the numb er density of these lines 
(see iDaU'Aglio et all [20081 12009D . With these synthetic 
Lya forest spectra, we illustrate both the optical depth 
decrease due to the LOSPE and the TPE of QSOs and 
then compare them with current observational results on 
type 1 QSOs in Section [H Considering that future X-ray 
observations potentially discover a large number of type 
2 QSOs, the TPE of type 2 QSOs is also investigated in 
Section [H Our simulation results demonstrate that the 
density enhancement near the QSOs, the QSO lifetime, 
and the opening angle of the tori associated with the 
QSOs can be simultaneously constrained through mea- 
surements on both the LOSPE and the TPE. Discussions 
and conclusions are given in Section [51 

In the paper. We adopt the Hubble constant Hq = 
71kms~^ an d the cosmological p arameters (r2m,^^A) = 
(0.27,0.73) (jKomatsu et al.|[2009f) . 



2. GEOMETRY OF THE PROXIMITY REGION OF A QSO 

In this section, we illustrate how the geometry of 
the proximity region of a fgQSO in the distant ob- 
server's frame is affected by the age of the QSO and 
the anisotropic feature of its radiation. In general, we 
describe the region near a fgQSO by the physical dis- 
tance of a point from the QSO R and its polar angle 6 
measured from the observer's LOS (i.e., OC in Figure [1]). 

From Lya forest spectra, the distant observer deduces 
the proximity region which can be affected by the fgQSO 
radiation through the light passing by the region. Con- 
sidering the finite speed of the light, the combination of 
the light traveling time from the affected region to the ob- 
server and the time of the fgQSO light traveling from its 
location to its affected region should not be longer than 
the light traveling time from the fgQSO to the observer 



by the QSO age. That is, in the distant observer's rest 
frame, photons emitted from the QSO can reach a phys- 
ical distance of at most R*{6^) = |0A| from the QSO 
as shown in Figure [1] and i?*(l — cos 6'*) — ctq, where 
Oif is the angle of OA from the observer's LOS OC, tq 
is the QSO age, and c is the speed of hght. The i?*(0*) 
gives the apparent time-delay envelope within which the 
ionization state of (hydrogen) atoms can be possibly af- 
fected by UV radiation from the QSO. 

Considering a bgQSO of which the LOS to the observer 
is separated from the fgQSO by a transverse proper dis- 
tance of i?^, the intersection point of the bgQSO LOS 
with the time-delay envelope (point A in Figure [ll is 
given by 6^ = 2 arctan(cTQ/i?_L). The Lya forest in 
the spectrum of the bgQSO may be then affected by the 
fgQSO at observational wavelength A < 1216A(1 -I- z^), 
where za is the redshift of the interaction point A. 
We have za — zig — (1 -I- Zfg)_B(2;fg)i?|| a/c, where Zfg 
is the redshift of the fgQSO, i?||,A = -t^^^ctq is 
the the physical distance from A to B shown in Fig- 
ure [1] point B is located on the LOS to the bgQSO and 
has the same distance/redshift to the observer as the 

fgQSO, and E{z) = Ho^/n^lT7fT^. The pho- 
toionization enhancement at a particular point on the 
LOS to the bgQSO due to the fgQSO can be charac- 
terized by the ratio of w = rgso/PuvB, where Fqso is 
the photoionization rate due to the QSO UV radiation 
and FuvB is the photoionization rate due to the cosmic 
UVB. The latest constraints on the cosmic UVB show 
that FuvB is roughly a constan t over redshift from 2 
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see also [Haardt k Madaul Il996l: iRauch et all 119971: 



McDonald fc Miralda-Escud^ 120011: iMeiksin fc Whitel 



20041: iTvtler et al.ll2004 IKirkman et al.ll2005[) . In this 
paper, we choose Fuvb — 10~-^^s~-^. The Fqso at 
a distance R from the fgQSO is ^^^^i^^^^TU^^u, 

where Lqso,i^(^) is the luminosity of the QSO per unit 
frequency at its intrinsic frequency v emitted at time 
t — Tq — R{1 — cos6)/c since its nuclear activity is trig- 
gered, 1^0 is the Lyman limit frequenc y, amii') j s 
the HI photoionization cross-section (|Osterbrocklll989D . 
and R can be approximated as the corresponding lumi- 
nosity distance as long as it is small enough. In this 
paper, we assume Lqsom oc (v/vq)^^-'^ . Due to the time 
delay, we have Fqso — and uj — at those points with 
R>R,i9). 

In the AGN unification model (e.g.. lUrrv &: Padovanil 
Il995f l. however, the UV radiation from a QSO can be 
highly anisotropic depending on whether the torus as- 
sociated with the QSO blocks its radiation along the 
LOS to the observer, and this geometrical effect results in 
the observational differences between type 1 and type 2 
QSOs. The anisotropic radiation of QSOs leads to differ- 
ent ionization states of its surrounding IGM in different 
directions and thus we may have different observational 
PEs for type 1 and type 2 QSOs as follows. 

• For a type 1 QSO, the UV radiation may be blocked 
in the transverse direction and the ionization state 
of the IGM in that direction may be not affected 
by the UV photons from the QSO. We denote the 
half opening angle of the dusty torus associated 
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Fig. 1. — Schematic diagram for the proximity region of a fgQSO that may be affected by the UV photons emitted from it. The fgQSO 
is at point O, and the distant observer is located in the direction OC. The ED is the direction from a background bright UV source or QSO 
to the observer. The left and right panels illustrate the cases for a type 1 fgQSO and a type 2 fgQSO, respectively. The intrinsic QSO UV 
radiation is assumed to be isotropic, but it may be shielded in some directions by a dusty torus as shown by the shaded area in both panels. 
The half opening angle of the torus is So = 45°, and the offset of the torus axis (the dot-dashed lines) from the direction OC is "i/o = 15° in 
the left panel and 105° in the right panel, respectively. The shaded region is rotationally symmetric around the torus axis. For illustration 
purpose, the LOS to the bgQSO is assumed to be specifically on the plane of the torus axis and the LOS to the fgQSO { see Appendix for 
discussion on more general cases); and it has a physical distance of i?x = from the fgQSO and interacts with the shaded region at 

points E and F. The thick solid curves represent the envelope surfaces of the proximity regions that can be affected by the fgQSO radiation 
for the QSO age of 3.2 X lO^yr, lO^yr, 3.2 X lO^yr, and 10''yr (from top to bottom), respectively. These surfaces are rotationally symmetric 
around the LOS to the fgQSO. The dashed circle represents the surface on which the photoionization rate enhancement due to the central 
QSO is w = rqso/ruVB = 1. Here the QSO is set to be at redshift ^fg = 2 with a Lyman limit luminosity of lO'^^'^ergs^^Hz"^, and the 
photoionization rate due to the UVB FuvB = lO^^^s^^. 



with the fgQSO by Oq (0° < Oq < 90°) and the 
torus axis offset from the LOS by an angle of 4'o 
(0° < "bo < 180°; see the left panel of Figure 
[T|). The intersections of the bgQSO LOS with the 
light cone confined by the torus associated with the 
fgQSO are E and F. Given a transverse distance 
R±, the proper distances from E to B and F to B 
are given by the two solutions (i?||,Ej ^||,f) of the 
equations governing the light cone and the cylin- 
der with a radius of R± (see details in Appendix), 
and correspondingly the redshifts of E and F are 
given by 2b ~ Zfg - (1 + ztg)i!;(zfg)i?||_E/c and zp ^ 
zig — + zig)E{z[g)R^\ p/c, respectively. Therefore, 
the UV photons from the fgQSO can only affect the 
spectrum of the bgQSO in the wavelength ranges of 
1216Amin(l-KzE, 1 + za) < A < 1216A(H-za) and 
A < 1216A min(l + ZA, 1 + ^f) as ze > zp (see the 
left panel of Figure [T|). Outside of these wavelength 
ranges, a; = 0. 

• For a type 2 QSO, the UV radiation is blocked from 
the view of a distant observer by a dusty torus. 
This QSO cannot be seen in the optical-UV band 
but is detectable in the hard X-ray. However, the 
UV radiation from a type 2 QSO may significantly 
alter the ionization state of the IGM in the trans- 
verse direction of its proximity region (see the right 
panel in Figure [1]). For a type 2 QSO, the UV pho- 
tons can affect the Lya forest spectrum of a bgQSO 



in the wavelength range of 1216A min(l + ze, 1 + 
za) < A < 1216A min(l + ZA,l + zp) (see the right 
panel of Figure [IJ. For those cases without so- 
lutions to the two equations, the spectrum of the 
bgQSO is not affected by the PE with uj = 0. 

We note here that the TPE of type 2 fgQSOs could 
be the same as that of type 1 fgQSOs if the obscura- 
tion, which leads to the classification of type 1 and type 
2 QSOs, is due to isotropically distributed clumpy ab- 
sorbers, rather than the uniform torus-like structure. 

3. MONTE-CARLO SIMULATIONS OF SYNTHETIC 
SPECTRA OF LYa FORESTS NEAR A FGQSO 

In principle, both the LOSPE and TPE can be 
well understood by detailed 3D numerical simula- 
tions that can succe ssfully reproduce the Lya forest 
at any redshift (e.g., [CroftI [20041: ISchirber et all 12004 
iFaucher-Giguere et all l2008bh ~ Those simulations with 
detailed radiative transfer may be able to accommodate 
all the statistical and systematic effects and the density 
enhancement in the vicinity of QSOs, but they are very 
time consuming and dependent on the assumptions of 
the QSO host dark matter halos. In this paper, we 
alternativ ely adopt a Monte Carlo m ethod, proposed 
by Worsc ck fc Wisotzkil (|2006l see also iDairAglio et al.l 
[20081 12009? ) ■ to generate Lya forest spectra for a large 
number of LOSs. This Monte-Carlo method has been 
demonstrated to be able to reproduce the Lya forest 
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spectra of QSOs similar to observed ones and be useful 
in measuring the UVB through the QSO proximity effect 
(jPairAglio et al.ll2008[ ). With detailed statistical con- 
sideration on different factors that may affect the prox- 
imity effect of the fgQSOs, including the QSO age, the 
anisotropic UV radiation and the density enhancement, 
etc. (as discussed in Section [5]), we simulate both the 
LOSPE and TPE for a large number of fgQSOs. 

The procedure adopted to generate the mock QSO Lya 
forest spectra is based on the observations that each Lya 
absorption line with Voigt profile can be described by the 
HI column density -/Vhi and the Doppler parameter b of 
its corresponding absorption gas and that the comoving 
number density of the absorption lines per unit z, -/Vhi, 
and b, denoted by n{z, N}ii,b)^ can be well described by 
the following three distributions: 

1. the redshift distribution, approximated by a power- 
law form with n oc (1 -|- z)'^ , where 7 = 



2.18 (iKim et al.l 120011: 


Schave et al.ll2003l 


iFaucher-Gieuere et al.l 


2008ch: 



2. the HI column density dis tribution with n oc iVjjj , 
where the power ^ ~ 1.5 (|Kim et al.ll200H l: 

3. the Doppler parameter distribution wit h n oc 
b~^ e xp[— 6^/6*^] where b^r — 24kms^^ (|Kim et al.l 
[2OOI . 

The column densities of the simulated absorbers 
are limited to be within the range lO^^cm"^ < 
A'hi < lO^^cm^^ and the Doppler parameters within 
lOkms""'^ < b < lOOkms"^. Those three distributions 
are assumed to be independent one another. For the 
purpose of this paper, we populate absorbers in each line 
of sight according to the above distributions and generate 
mock samples of the Lya forest spectra. Gaussian noise 
is added to the simulated Lya forest spectra and the 
spectral resolution is also considered in order to match 
the signal-to-noise level of those observed samples dis- 
cussed in Section 4. For simplicity, the emission lines 
are not added to the spectra of mock QSOs. Note that 
the adopted Monte-Carlo method naturally includes the 
Poisson variance of the nu mber of absorption line s in the 
near zone of fgQSOs (e.g.. iDall'Aglio et al.l[2008[ ). 

In the QSO near zones, the mass density may be signif- 
icantly enhanced relative to the cosmic average density, 
i.e., (A(i?)) = {p{R)) Ip > 1, where (A(i?)} and {p{R)) 
is the mass overdensity and the mass density at a dis- 
tance R from central QSOs averaged over a large sample 
of QSOs with similar properties, respectively, and p is 
the cosmic average density. A simple photoionization 
equilibrium model can give the neutral hydrogen den- 
sity of an absorber nni oc 5^~°-'^^, where S is the mass 
overdensity of the absorber S and rj ^ 0.62 is the in- 
dex of the power-law temperatu re-density relation for 
the Iqw-temperature IGM (e.g.. iFaucher-Giguere et al.l 
I2008hl : iSchavd [2OOII: iHui fc GnedinI 1199711 . Assuming 
that the column density distribution of Lya absorbers 
in the QSO near zones (with mean density enhancement 
factor (A(i?))) also follows a power law oc N^f with 
the same slope /3 as that for the cosmic average (item 2 
listed above) , the column density of each absorber gen- 
erated from the above distribution (item 2) should then 
be replaced by A^hi (A(i?))^^"'^'' / (1 +0;), considering of 



both the effects of the enhanced density and UV ion- 
izing flux in QSO near zone s (cf., see Equation 30 in 
iFaucher-Giguere et al1l2008b[ ). Note that the effect of 
that replacement is equivalent to the effect of increasing 
the number of absorbers proportionally. 

As illustrated in Figure [1] different choices of the 
fgQSO age, the half opening angle Oq, and the offset 
angle of the torus associated with the fgQSO will 
affect the proximity region that the fgQSO can illumi- 
nate. In our simulations, we set seven different values 
for the lifetime of QSOs Tit, and the QSO age tq is 
either fixed or randomly chosen over the range (0, rit). 
We choose six sets of values for the opening angle of 
the torus, i.e., Oq = 0°, 30°, 45°, 60°, 72°, and 89°. The 
6o = 0° corresponds to the cases of no ionization PE 
from the central QSOs, and the ©o — 89° corresponds 
closely to no obscuration to the UV radiation from the 
fgQSOs. The other specific values of ©0, i.e., 30°, 45°, 
60°, and 72°, are chosen as their corresponding abun- 
dance ratios of type 2 to type 1 QSOs (6.5, 2.4, 1, and 0.4, 
respectively) ar e roughly in the range estimated by ob- 
servations Ce.g.. iTreister fc Urrvl[200g: iGilh et aH 120071: 
iTreister et al.ll2010l ). We also randomly choose within 
the range from ©o to 180° — ©0 for type 2 QSOs, and 
other range for type 1 QSOs. 

According to the above settings, we first simulate the 
Lya forest spectra of QSOs affected by its own UV radia- 
tion to account for the LOSPE. In order to compare with 
the obse rvational results on the L OSPE (and the TPE 
later) bv lKirkman fc Tvtleil (12008) (hereafter KT08) ob- 
tained from a sample of 130 QSO pairs, the observational 
redshift of the fgQSOs and its luminosities at the Lyman 
limit frequency in our simulation are set to be Zfg = 2 
and ijyp = 5 X 10^°erg s"^ Hz^^, which are roughly the 
mean redshift and luminosity of the sample in KT08. 
Note that most of the QSOs in KT08 sample are in the 
redshift range z ^ 1.8 — 2.6, and the redshift dependence 
of the optical depth has been corrected in KT08. As the 
actual halos hosting these fgQSOs and the overdensity 
distribution surrounding them are not well known, we 
only consider the average effects of the sample in this 
paper, not going to the detailed redshift and luminosity 
distributions of QSOs as shown in the KT08 sample. We 
first create 100 realizations of 130 independent synthetic 
spectra to check whether the sample variance could be 
significant in the interpretation of the observational re- 
sults. And then we also create larger mock samples with 
500 independent synthetic Lya forest spectra to check 
whether the QSO properties can be extracted effectively 
from a sample with 500 spectra or more. Similar to 
KT08, the amount of absorption in each spectrum of the 
mock samples is quantified by DA = 1 — F/C, where F is 
the flux and C is the continuum level, and a uniform flux 
decrement due to the metal absorption DA^etai — 0.025 
is also added to each spectrum at wavelength higher than 
the Lya emission line. We also simulate the Lya forest 
spectra of bgQSOs, which are affected by the UV radi- 
ation from fgQSOs, to account for the TPE. The trans- 
verse distances of fgQSOs from the LOS to their bgQSOs 
R± are randomly chosen from 0.1 Mpc to 3 Mpc, within 
the range of the observational sample in KT08. We also 
create 500 independent synthetic spectra for each set of 
parameters (rage, ©0) for both type 1 and type 2 fgQSOs, 
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respectively. For each mock sample, we then obtain the 
average flux decrement (DA) of the 500 spectra. Com- 
parison of the simulated LOSPE and TPE with the ob- 
servations may thus put some constraints on the intrinsic 
properties of QSOs. 

Similarly, we create 500, if not otherwise specified, in- 
dependent synthetic spectra of a QSO/bgQSO, which are 
affected by either its own UV radiation or the UV ra- 
diation from a fgQSO at Zfg = 4 with a mean Lyman 
limit luminosity of L^^ = 8 x lO'^^erg s^^ Hz~^, for each 
set of parameters (rn, Qq)- The LOSPE obtained from 
the set of synthetic Lya forest spectra can be used to 
compare with that obtained by iGuimaraes et al.l ()2007( ) 
for a sample of QSOs with simil ar redshift and mean Ly - 
man limit luminosity. Similar to ^Guimaraes et al.l (|2007t ). 
the absorption here is quantified by the median of the 
optical depths r = - ln(F/C) of aU the LOSs. With 
so me constraints obtained from the LOSPE observations 
bv IGuimaraes et"al1 ()2007| ). we also obtain the expected 
TPE of type 1 and type 2 fgQSOs from the mock sam- 
ples. 

4. SIMULATION RESULTS 

We illustrate the simulation results on both the 
LOSPE and TPE for both the samples with relatively 
low luminosities at low redshift (z ^ 2 and L^^ ~ 
5 X 10^°erg s^^ Hz^^; Section |4T|) and the samples with 
high luminosities at high redshift (z ~ 4 and L^^ ~ 
8 X lO'^^erg s~^ Hz~^; Section . In order to compare 
with those observations by KT08 and IGuimaraes et al.l 
()2007( ). the spectra are averaged over resolution of 
and O.SA, and the signal to noise ratio are set to be 
11.2 and 25 for the low-luminosity and high-luminosity 
sampl es similar to those in KT08 and Guiniaraes et al] 
()2007| ). respectively. Along with the detailed presenta- 
tion on how the detected PEs can be affected by the 
QSO lifetime, the geometric structure of the torus, and 
the QSO environment below, the effect of different QSO 
luminosities is also revealed. 

4.1. Low-luminosity samples at z ^ 2 

4.1.1. The LOSPE of type 1 fgQSOs 

Figure [2] shows the average absorption near 
QSOs/fgQSOs along their LOS with the origin corre- 
sponding to their redshift. The histogram with error 
bars shows the DA measured by KT08. The color lines 
and symbols in panels (a) and (b) show the results 
obtained from 100 realizations of a mock sample with 
130 synthetic Lya forest spectra and the results from a 
mock sample with 500 synthetic spectra, respectively. 
The value plotted for each bin is the mean value of DA 
from the Lya forest spectra, and the error bar is taken to 
be standard deviation of the 100 realizations in panel (a) 
and the standard error of the mean in panel (b), respec- 
tively. The color lines and symbols show the expected 
DA estimated by assuming the following several cases: 

(1) on average the density in the QSO near zones is not 
enhanced (blue line and points), i.e., (A(i?)) — 1; or 

(2) there is no UV radiation from the central QSOs and 
on average the density near the QSO is not enhanced, 
i.e., uj — {) and (A(i?)) = 1 (green line and points); 
or (3) there is a density enhancement near the QSO 
described by (A(i?)) ~ 1 + Co(i?/C/?)-P exp[-(i?/Ci?)«] 



with Co = 2.2, Cr = S.OMpc, p = 0.7, and q = 2.0 
(red line and points)^. As shown in panels (a) and (b), 
the histogram is systematically and substantially higher 
than the blue line but consistent with the green line 
(no LOSPE), which suggests that the sample variance 
is highly unlikely to be responsible for the excess of 
absorption in the QSO near zone and the density in the 
QSO near zones is enhanced. Indeed, the expected DA 
can be consistent with the observations if the density is 
enhanced in the QSO near zone (as shown by the red 
line and points in panels a and b). Note also that the 
error to the mean of DA for a sample of 500 spectra 
is sufficiently small, which suggests that the physical 
properties of QSOs can be extracted from the PE for 
samples with 500 spectra or more. 

The PE analysis is sensitive to the accuracy in the 
systemic redshift estimation of each QSO/fgQSO being 
studied if the size of its proximity region is smaller than 
or comparable to the scale corresponding to the systemic 
redshift error. However, the systemic redshift estimated 
from the UV emission lines of a QSO is generally sys- 
tematically smaller than the systemic redshift of its host 
galax y by more than a few hundred kilometer per sec- 
ond. iRichards etHI (|2002D find that CIV gives a red- 
shift smaller than Mg II by 824 kms~^ with a scat- 
ter of 511kins~^ fr om a subset of SDSS QSO spectra. 
iNestor et ahl (|2008D find that the [OIII], which should 
better represent the systemic redshift, gives a systemic 
redshift larger than Mg II by 102 kms~^. According 
to these results, the errors in the systemic redshift es- 
timates of QSOs can be divided into two parts, i.e., the 
mean systematic offset and the scatter around the mean 
offset, which are denoted as the "systematic error" and 
the "random error" belo w, respectively. Ac cording to 
IRichards et"all (|2002?) and 'N estor et al.l(|2008D . we adopt 
a systematic error of 926kms~^ in the estimates of the 
fgQSO systemic redshifts and a random error with dis- 
persion 511kms~^ (as the redshift of the majority of the 
QSOs in KT08 is estimated by CIV), unless otherwise 
specified. Note that the systematic error adopted in 
KT08 is smaller, ~ 753kms~^, and the dispersion of the 
random errors is 676kms~^, and we also try those er- 
rors in our simulations below (i.e., in Figure El where 
the systemic redshift of each mock QSO is set to be a 
value of Zfg = 2 plus a Gaussian distributed random er- 
ror). In order to compare with the observational results 
from KT08, an additional correction in the systematic 
error of (^Zoffsct) — 173kms~^ needs to be taken into ac- 
count, i.e., the systemic redshift of each QSO host in the 
mock samples is set to Zfg -I- ^ZoSsct, where Zfg is set to 
2 and the random error in 5zoSsct is also assumed to be 
Gaussian distributed (e.g., in Figures |3Jl, [H [SI 13 be- 
low). For comparison, a zero error of ^ZoSsct is assumed 

^ Note that this assumed form of (A(_R)) follows the tendency 
of approaching 1 at sufficiently distant regions and increasing with 
the decreased R. Here (A(i?)) can be taken as the effective density 
enhancement factor but not the exact physical overdensity of each 
absorber at R. The physical overdensity of an absorber at R is 
(5 oc (A(_R)). The detailed values used for the involved parameters 
come empirically from the comparison of our simulation results 
with the observations below. While this should be sufficient for the 
purpose of this paper, those constraints may be improved by future 
high-resolution spectra of QSOs and the PEs to be revealed and the 
comparison with cosmological simulations of PEs with including 
detailed 3D radiative transfer processes. 
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Fig. 2. — The line of sight proximity effect. The histogram 
in panel a (or b) shows the measured flux decrement DA at a 
proper distance R|| along the LOS from the central QSO, averaged 
over an observational sample of 130 QSOs with redshift ~ 2 and 
Lyman limit luminosity L„g ~ 5 X lO^^erg s~^ Hz~^ (adopted 
from KT08). The color lines and points in panel (a) represent 
the mean DA obtained from 100 realizations of a mock sample 
with 130 synthetic Lyo forest spectra, where the error bar for each 
bin is the standard deviation of the sample, while in panel (b) they 
represent that obtained from mock samples with 500 synthetic Lyo 
forest spectra, where the error bar is the standard error on the 
mean. The blue line and solid circles show the simulated DA due 
to the PE of the central QSOs if the average density in the QSO 
near zones is similar to the cosmic average; the green line and 
solid circles show the simulated DA if there is no proximity effect 
due to the UV radiation from the central QSO and no density 
enhancement in the QSO near zones; and the red line and solid 
circles show the simulated DA if the density is enhanced in the 
QSO near zones effectively by a factor of (A(iJ)} (see details in 
Section l4.1.1l l. All the circle points are centered on the same bins 
of as the observational data but have been slightly offset to the 
right for graphical clarity. 

in Figure |3^-c below. Generally, the larger the system- 
atic error in the systemic redshift estimates, the less the 
enhancement of the density in the fgQSO near zones that 
is required to reproduce the observed LOSPE (see Fig- 
ure [2]). 

4.L2. The TPE of type 1 fgQSOs 

Assuming that the effective density enhancement in the 
near zones of the fgQSOs is the same as the one required 
to reproduce the observations on the LOSPE by KT08^, 
we generate mock samples with 500 synthetic Lya for- 
est spectra and obtain the TPE effect from these sam- 
ples. Figure [3] shows the expected TPE obtained from 
the mock samples with different settings of the QSO age 
tq, the opening angle of the torus Oq, and the error 
in the estimation of the fgQSO systemic redshift. Fig- 

^ We assume that the density enhancement in the QSO near 
zone is isotropic. We do not consider an anisotropic distribution 
of the density enhancement in this paper. 



ure shows that for a small QSO age tq, the excess 
of DA is significant near the fgQSOs because of the sig- 
nificant density enhancement in the immediate vicinity 
of the fgQSOs. This excess of DA decreases with in- 
creasing Tq because the significance of the suppression 
of the absorption by the UV photons from the fgQSOs 
increases with increasing tq. In principle, the depen- 
dence of the TPE on the QSO age tq, as shown in Fig- 
ure |3K, suggests that the TPE can be used to constrain 
the QSO lifetime once the density enhancement in the 
fgQSO near zones is determined by the LOSPE. The dif- 
ferences among large tq cases shown in the panel are 
small because the time during which photons cross the 
proximity regions are smaller than or at most compara- 
ble to the QSO age in these cases. However, none of the 
simulated DA (even the one with the shortest tq) in Fig- 
ure can match the observational results on the TPE 
obtained by KT08. 

Figure [Sb shows the dependence of the TPE on the 
half opening angle Oq of the tori associated with the 
fgQSOs. For the extreme case of Oq = 0° (without the 
PE due to the fgQSO UV radiation), the excess of the 
DA due to the density enhancement near fgQSOs is the 
most significant (red line and points). With increasing 
8o, the excess of DA becomes less and less significant 
because the region that can be affected by the UV pho- 
tons escaping out from the central engine becomes larger 
as shown by the color lines and points. For 9o — 89° 
(almost without obscuration to the fgQSO photons; ma- 
genta line and points), the expected excess of DA due to 
the density enhancement is balanced significantly by the 
proximity effect due to the fgQSO UV radiation. 

Figure [3t shows the dependence of the TPE on the 
lifetime of fgQSOs, in which the half opening angle of 
the torus is fixed to 6o = 60°. This value of 6o is 
roughly in the range d etermined by obs e rvations (e.g., 
Treister fc Um^ 120061 : iGilh et al.l l20Q7t iTreister et atl 



2010f ). In this panel, the value of the QSO age tq is cho- 



sen randomly over a range from to the QSO lifetime 
Tit, instead of being a constant as labeled for each line in 
panel (a). Because of the obscuration in the transverse 
direction to the UV radiation from the fgQSOs, there is 
significant DA excess even for the case of T\t = 10^'^ yr. 
Compared to the case without obscuration shown in Fig- 
ure [3^, the dependence of the DA excess on th becomes 
less obvious if the obscuration to the UV radiation from 
the fgQSOs is significant. 

Note that one characteristic timescale for the lumi- 
nosity evolution of a QSO is the Salpeter timescale Tgp, 
which is ~ 10^ '' yr if the mass-to-energy conversion ef- 
ficiency of its nuclear activity is ~ 0.1 and the Edding- 
ton ratio is 1. If the age of a fgQSO is larger than the 
Salpeter timescale, the proximity effect may be substan- 
tially less significant for regions faraway from the fgQSOs 
with i?(l — COS0,) > cTsp, considering the luminosity evo- 
lution of the fgQSOs may be significant. 

Figure [3ji shows the expected TPE obtained from the 
mock samples after taking into account the errors in 
the systemic redshift of the fgQSOs. As in panel (c), 
the opening angle of the tori associated with the fgQ- 
SOs is fixed to 60° and the age of fgQSOs tq is ran- 
domly chosen over a range of (0, Tit). As seen from 
the panel, the observations by KT08 can be well rcpro- 
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Fig. 3. — The transverse proximity effect for type 1 foreground QSOs at redsiiift ~ 2. Tiie histogram shown in each panel represents 
the average flux decrement DA obtained by KT08 for a sample of background QSOs as a function of proper distances along the sightline, 
with the origin corresponding to the redshift of the paired fgQSOs; and other lines and symbols represent our simulation results. In panel 
(a), the simulation results of DA is obtained under the assumption that there is no obscuration to the fgQSO UV radiation and no error in 
the estimates of the systemic redshifts of fgQSOs. Different colors represent different QSO ages, as labeled in the panel. Panel (b) shows 
the simulation results obtained by assuming that the QSO lifetime is lO'^'^yr, and the QSO age is randomly chosen over a range from to 
each assumed QSO lifetime. Different colors represent different half-opening angles of the torus associated with each fgQSO Sq. Panel (c) 
shows the simulation results obtained with ©o = 60° and different values of QSO lifetime. All the simulations in panels (a)-(c) assume a 
zero error in the estimates of the fgQSO systemic redshift. Panel (d) show the simulated DA obtained by the same assumptions as those 
in panel (c), except for accounting for both the systematic and random errors in the estimates of the fgQSO systemic redshift. All the 
circle points shown in the figure are centered on the same bins of as the observational data but have been slightly offset to the left for 
graphical clarity (similarly in Figures [5] and [6] below) . 



duced if the lifetime of fgQSOs is ^ 10^-^ yr (green line 
and points). Our calculations also show that the ex- 
pected DA for cases with ru < 10^ yr cannot match the 
observations for any given Qq, which is different from 
the suggestion that the fgQSOs have had their current 
UV luminosities for less than approximately a million 
years made in KT08. Our results suggest that the life- 
time of the fgQSOs is consistent with being a few times 
10'' yr and a substantially smaller lifetime (< 10^ yr) 
is not required by the PE, and it is also consistent 
with those constraints obtained by the detections of the 
TPE through He II Lyct for e st or metal absorptio n lines 
(e.g.. IJakobsen et al l [200l iWorseck fc Wisotzkil [200l 
IWorseck et al.ll2007t iGon^alves et al.ll2008D 



Note here that the two components of some QSO pairs 
in KT08 are at similar redshift and therefore the TPE 
of those fgQSOs may be affected by the bgQSOs. The 
UV light from the bgQSO near a fgQSO may contribute 
some to the ionization of the near zone of the fgQSO, 
but this effect should be more significant at the back- 
side than the front side of the fgQSO and lead to a more 

significant decrease in the DA at ^ —3 5 Mpc. 

This effect cannot explain the significant excess of DA 
at ^ —3 - —5 Mpc but rather require a more sig- 
nificantly enhanced density in the fgQSO near zone. In 
addition, if the density is also enhanced in the near zone 
of the bgQSOs, a more significant excess of the absorp- 
tion at Rh ^ —5 Mpc can be expected, which might bring 
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the observed DA at R\\ — —5 Mpc in better consistency 
and strengthen the above conclusions. 

In order to check the significance of the sample vari- 
ance on the TPE obtained by KT08, we again generate 
100 realizations of mock samples, each with 130 synthetic 
Lya spectra. Figure |4] shows the expected TPE for those 
mock samples with similar settings to that shown in Fig- 
ure [3ji. According to this Figure, the KT08 results are 
still consistent with that the QSO lifetime is as long as 
a few 10^ yr and the torus half opening angle Qq ^ 60°. 
For the cases of T\t < 10® yr, the probability that the 
inconsistency between the KT08 TPE results and the 
expected TPE are simply due to sample variance is low. 

There are still some uncertainties in the current esti- 
mation of UVB which may affect the simulation results. 
In our calculations, we set the UVB as Fuvb = lO^^^s^^, 
w hich might be somewhat large as the latest estimation 
of iFaucher-Giguere et afl (|2008aD is only 0.5 x lO'^^g-i 
If the UVB is set to be this smaller value, a more sig- 
nificant LOSPE would be expected, and the effective 
density enhancement factor {A{R)) should be slightly 
larger compared to that in Section 14.11 in order to fit the 
LOSPE estimated by KT08. Therefore, the QSO life- 
time Tit and/or the torus opening angle 8o are required 
to be even larger than that given in Section [4.1l in order 
to reproduce the TPE. But if the UVB is unreasonably 
much larger than 10~^^s~^, the excess of absorption at 
the transverse direction detected by KT08 cannot be ex- 
plained by simply changing r^, 6o and (A(i?)). 

We caution here that other uncertainties could also af- 
fect the results presented here quantitatively. First, the 
DA estimated from observations may be affected by the 
continuum fitting. However, the error in the continuum 
fitting is typically on the percentage level de pending on 
the signal-to-noise ratio of the QSO spectrum ()Kim et al.l 
[2007). which is not likely to change the observational re- 
sults by KT08 on PE qualitatively. Second, the combined 
sample in KT08 are obtained from several different in- 
struments and is highly heterogeneous. The ignoring of 
the detailed exact redshift and luminosity distributions 
of the QSO pairs in our simulations is sufficient for the 
demonstration purpose in this paper, but it may intro- 
duce some uncertainties to the resulted PE. Future works 
on extracting QSO properties from the PEs should con- 
sider the uncertainties. In addition, the luminosity of 
QSOs may evolve or fluctuate on timescales 10® yr, which 
may lead to more sign ificant excess of a bsorption in the 
transverse directions (| Adelberged l2004 ) . However, such 
a luminosity variation of QSOs is not required according 
to our simulations. 

Figure [5] shows the expected DA obtained from the 
mock samples by assuming that both the systematic and 
random errors in the estimates of the fgQSO systemic 
redshifts through CIV are the same as those adopted in 
KT08. Compared to the density enhancement used in 
Figure [2l here a larger value of it is used so that the 
LOSPE of fgQSOs obtained by KT08 can also be re- 
produced well (see the red line and points in Figure (5^). 
With the density enhancement required by the LOSPE, 
the DA distribution near the fgQSOs indicated by the 
bgQSO spectra (i.e., the TPE) is also calculated as shown 
in Figure [5}d. We find that the observational asymmetric 
distribution of the excess in DA near ^ 0, i.e., the 
lack of excess in DA at R\\ ^ 2.5 Mpc and the significant 
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Fig. 4. — The transverse proximity effect of fgQSOs expected 
from 100 realizations of mocls samples each with 130 synthetic Lyo 
forest spectra. Legends are similar to that for Figure [3jl, except 
the errorbar to each bin is the standard deviation of the sample. 
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Fig. 5. — The line of sight proximity effect and the transverse 
proximity effect of fgQSOs that both the systematic and random 
errors in their systemic redshift estimation used are the same as 
that in KT08. Legends of panel (a) are similar to those of FigurelSl 
and legends of panel (b) are similar to those of Figure [Sjl. 

excess of DA at ^ —5 Mpc, cannot be simultaneously 
re-produced for any given rn and Oq, in contrast to the 
results in Figures |3Ji and|31 

4.1.3. The TPE of type 2 fgQSOs 

Assuming that the effect of the density enhancement 
near type 2 fgQSOs is the same as that of type 1 fgQ- 
SOs with similar intrinsic L^^ constrained by the LOSPE 
(see Figure [2]), we generate synthetic Lya forest spectra 
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of type 1 bgQSOs whose light passed by the proximity 
region of type 2 fgQSOs to study the TPE of type 2 fgQ- 
SOs. Figure m shows our simulation results on their DA 
distribution near the fgQSOs. Figure [6^ shows the de- 
pendence of the TPE on the half opening angle of the 
torus 6o for the mock samples of type 2 fgQSOs, given 
the QSO lifetime. As seen from the panel, generally the 
larger the Oq, the less the excess of the absorption near 
^ Mpc. Figure [BJd shows the dependence of the 
TPE on the QSO lifetime, given the half opening an- 
gle ©0 — 60°. As seen from the panel, the differences 
in DA are relatively small among the cases with large 
(- 10^ - 10^-^ yr). If Tit ~ 10^ yr, the absorption 



Tit 



at ^ Mpc is larger than those cases with large rit 
due to the time-delay effect. Generally the absorption 
in the region closer to the observer > 0) should be 
relatively large compared to that for type 1 fgQSOs (see 
Figure [Sji) , because the near side is more likely to be 
obscured from the fgQSOs. That effect, though weak, 
is shown in Figure [7] (see some relative deep dips of the 
curves at R\\ ~ 2.5 Mpc). This contrast between the 
absorption DA of type 2 fgQSOs and that of type 1 fgQ- 
SOs appears more significant for higher luminosity mock 
samples, as will be seen in Section 14.2.31 and Figure [Til 
below. 

Note that in this paper, the UV radiation is assumed to 
be intrinsically isotropic, after removing the anisotropic 
effect due to the torus. However, if the UV radiation 
is intrinsically anisotropic, the TPE of type 2 fgQSOs 
may be different from those expected from the above 
calculations. For example, if the UV (and X-ray) radia- 
tion is relatively stronger along the direction of the torus 
axis than in the direction along the torus plane, the ex- 
pected DA excess for type 2 fgQSOs should be smaller 
at ~ Mpc than those calculated above. 



4.2. High-luminosity samples at z ^ A 
4.2.1. The LOSPE of type 1 fgQSOs 

Figure [8] shows the optical depth distribution near 
QSOs/fgQSOs along their LOSs with the origin corre- 
sponding to their redshift (similar to Figure [2]). The 
red points with e rrorba rs show the results obtained by 
iGuimaraes et al.l ()2007D for an observational sample of 
^ 50 high-redshift QSOs [z ^ 4) with mean Lyman 
limit luminosity of L^^ ~ 8 x lO^^erg s~^ Hz~^. The 
blue open circles represent the median optical depth ob- 
tained from our mock sample with 500 spectra affected by 
fgQSOs by assuming that the density distribution near 
the fgQSOs is the same as that of the cosmic average, 
which are clearly offset from the observations. Note that 
the errors in the estimates of the fgQSO systemic red- 
shift are also considered in our simulati on here by adopt- 
ing th e method same as that used in IGuimaraes et al.l 
(|2007f) . i.e., the error is randomly chosen over the range 
of — 1500kms~^. The decrease of the simulated opti- 
cal depth with decreasing distance to the QSOs appears 
more significant than the observation results, as the ef- 
fect of the density enhancement in the proximity regions 
of these QSOs is ignored. Assuming that the density 
in the QSO proximity regions is enhanced by a factor 
of (A(i?)) ~ 1 + Co{R/Cr)'P ejip[-iR/CR)'J], the ob- 
servations can be reproduced by simulations (blue solid 
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Fig. 6. 



The transverse proximity effect for type 2 fgQSOs 



with the same intrinsic luminosity L^^ at redshift ~ 2 as those 
in Figure [3^3 and d. Legends for panel (a) and (b) are similar to 
those of Figure [Sji and d, respectively. As a reference, the excess 
absorption near type 1 fgQSOs measured by KT08 is also shown 
by the histogram here. 
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Fig. 7. — The difference of the transverse proximity effect be- 
tween type 1 fgQSOs and type 2 fgQSOs shown in Figure [3jl and 
Figure [B]]. The flux decrement in the proximity region of type 1 
fgQSOs and type 2 fgQSOs are represented by DA{T) and DA{IT), 
respectively. Legends are similar to those in Figure [3jl. 



circles) if Co = 0.9, Cr = 11 Mpc, p = 0.7, and q = 0.5.^ 
For this case, we set the number of spectra in each real- 
ization to be 50 and we simulate 100 realizations. Thus 

^ With the adopted parameters here, the effective density en- 
hancement is almost the same as that required for the low- 
luminosity sample of KT 08 at R < 3 Mpc but declines slower 
at J? > 3 Mpc (see Section |4]TJ . 
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Fig. 8. — The line of sight proximity effect for typo 1 QSOs. 
The red points show the optic al depth distribution in the QSO 
proximity region obtained by G uimaraes et aLl l |2007l ') for an obser- 
vational sample of QSOs with a mean Lyman limit luminosity of 
'~ 8 X lO'^^erg s~^ Hz~^ at redshift ~ 4. The blue open circles rep- 
resent the optical depth distribution obtained from our simulated 
mock sample by assuming that the density in the QSO proxim- 
ity regions is the same as that of the cosmic average. The blue 
solid circles with errorbars represent the simulated results by as- 
suming that the density distribution in the QSO proxi mity r egion 
is enhanced by a factor of A(iJ) (see details in Section l4. 2.11 1. 

the mean median optical depth and its standard devia- 
tion can be estimated from these realizations. As seen 
from Figure [HI the observations can be well matched by 
t he simulation resu l ts (bl ue solid points). 

iGuimaraes et alj ()2007t ) also suggest that the density 
enhancement in the near zones of QSOs is required in or- 
der to explain their measurements on the LOSPE. Com- 
pared with their estimates on the density enhancement 
in the near zones of those high luminosity QSOs, our 
estimates are smaller by a factor of 1.5 — 5 at a dis- 
tance of < 10 Mpc-1 Mpc but similar at a distance 
of > 15 Mpc. The main reason for this difference is 
that the Poisson variance of absorption lines in different 
sight lines to the near zones of QSOs, which is included 
in the Monte-Carlo simulations here, leads to an aver- 
aged absorption that is larger than a simple estimation 
obtained without consi dering of the Poisson variance (see 
IGuimaraes et al.ll2007[ ). Therefore, the density enhance- 
ment requir ed here is less significant compared with that 
obtained bv IGuimaraes et'al ? (2007) (see also discussions 
on the Poisson variance in IbairAglio et al. 2008,). 

4.2.2. The TPE of type 1 fgQSOs 

Although so far there is no observational mea- 
surement on the TPE of QSOs (with luminosity ~ 
8 X lO'^^erg s~ ^ Hz~ ^) similar to the sample in 
IGuimaraes et al.l ()2007[ ). we illustrate here the effects 
of different QSO properties on the TPE and demon- 
strate that both Tit and Oq can be simultaneously con- 
strained by the TPE if the density enhancement in the 
fg'QSO proximity region has been constrained by the 
LOSPE. We first assume that the effect of the density 
enhancement in the proximity regions of those QSOs is 
the same as that indica ted by the LOSPE obtained by 
IGuimaraes et al.l (|2007[ ) (see Figure [5]). Under this as- 
sumption, we synthesize a large number of Lya forest 
spectra of bgQSOs with different choices of the QSO life- 



time and the half opening angle of the associated torus 
(see parameter settings in Section [5]) . We then extract 
the median optical depth distribution in the fgQSO prox- 
imity regions from these spectra and show the results in 
Figure |H1 As seen from Figure [H^, given the half open- 
ing angle of the torus associated with the fgQSO, e.g., 
80 = 60°, generally the optical depth in the proximity 
region to the QSO (e.g., < 10 Mpc) decrease with 
increasing the QSO lifetime; and the decrease is signifi- 
cant initially at the frontside to the QSO, and becomes 
significant at its backside (e.g., ^ —4 Mpc) when 
the QSO lifetime th is long enough so that the region 
can be reached by the radiation from the fgQSO. That 
tendency is manifest especially for large Oq, as (1) for 
the case with short QSO lifetime (e.g., rit 10^ yr; 
see Figure injs), the insignificant decrease of the optical 
depth is not sensitive to different choices of &o; and (2) 
as shown in Figure |9]; and d, for large rjt, the decrease of 
the optical depth at < Mpc increases with increas- 
ing large Oo. As shown in the figure, the QSO lifetime 
and the half opening angle of the torus affect the opti- 
cal depth curves in different ways, so that they can be 
simultaneously constrained by the LOSPE and TPE of 
QSOs if Tit and 9o are not too small (e.g., rit > 10^ yr 
and 80 30°). Note here that the standard deviation, 
on the order similar to that shown in Figure O is small 
enough so that the constraints on the QSO lifetime and 
the torus half opening angle may be accurately extracted 
through the PEs of a sample with a few hundreds of high 
luminosity QSOs. 

Considering the possibility of inaccuracy in the errors 
of the QSO systemic redshifts used above, we also test 
that our results obtained above remain the same qual- 
itatively, if the systemic redshift errors is set to be the 
same as that for the low-luminosity low-redshift sample 
in Section 14.1.21 The overall shape of the optical depth 
distribution around ~ Mpc does not change signif- 
icantly, though it may shift slightly toward the negative 
direction, as discussed in Section B. 1.21 

4.2.3. The TPE of type 2 fgQSOs 

Assuming that the effect of the density enhancement 
near type 2 fgQSOs is the same as that of type 1 fgQSOs 
with similar L^^ ^ 8 x lO'^^erg Ilz~^ at z ^ 4, we 
show our simulated TPE of these type 2 fgQSOs in Fig- 
ure [ini As seen from Figure [TUk . given the half opening 
angle of the torus (e.g., 80 = 60°), the larger the QSO 
lifetime Tit, the smaller the optical depth at ~ Mpc 
because the UV photons from the fgQSOs are more likely 
to propagate to more distance regions in the transverse 
direction. Given the QSO lifetime, the larger the half 
opening angle of the torus, the smaller the optical depth 
at > Mpc (Figure [11)). 

Figure [11] shows the difference between the optical 
depth distribution near type 2 fgQSOs and that near 
type 1 fgQSOs, assuming that other parameters of these 
two types of fgQSOs (e.g.. Tit and 80) are the same. 
In contrast to the cases considered for fgQSOs with 
ij^o ~ 5 X 10^°erg s~^ Hz^^ at redshift '-^ 2, the differ- 
ence here is obvious at the near side of the fgQSOs and 
for some cases in the far side because the fgQSOs in the 
mock samples here are much more luminous than that in 
Section l4.1.3l Given the half opening angle of the torus, 
80 = 60°, the optical depth difference at ^ 0—15 Mpc 
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Fig. 9. — The expected transverse proximity effect of type 1 fgQSOs at redshift ~ 4 with Lyman hrnit luminosity of Lug ^ 8 X 
10'^^ erg Hz~^. Panel (a) shows the expected optical depth distribution in the proximity region of fgQSOs, in which the half opening 
angle of each fgQSO is assumed to be 60° and different curves represent the results obtained for different sets of the fgQSO lifetime 
Tit. The other three panels show the dependence of the optical depth on the half opening angle of the torus, given the fgQSO lifetime 
nt = 10^-^, 10^, 10^ '^, respectively 

increases with increasing tu but saturates when th > of type 1 (and/or type 2) fgQSOs. As the fgQSOs in 
a few times 10^ yr; while the difference at the far side these samples are much more luminous, their TPE is ob- 
< only becomes visible when T\t > 10^ yr. Given vious, without being smeared out by the errors in the 



a QSO lifetime, the optical depth difference is close to 
for either small half opening angle of torus Oq — >■ 0° or 
large Oq — > 90°; however, this difference is the largest for 
a medium Qq ~ 60°. That behavior of the optical depth 
difference is the combined effects due to the geometrical 
nature of the torus and the time delay of photons that 
propagated to the backside of fgQSOs. The sharp con- 
trast of the TPE due to type 1 fgQSOs from that due to 
type 2 fgQSOs should be useful to probe the properties 
of QSOs through the PE. 

The simulated results illustrated ab ove for the fgQ- 
SOs s imilar to that in the sample of iGuimaraes et all 
()2007[) again suggest that the QSO properties, such as 
the density enhancement in its proximity region, the life- 
time and the opening angle of the associated torus, can 
be constrained simultaneously by the LOSPE and TPE 



systemic redshift estimates of the fgQSOs. 

5. CONCLUSIONS AND DISCUSSIONS 

In this paper, we have investigated both the LOSPE 
and the TPE due to type 1 QSOs and the TPE due to 
type 2 QSOs. To do so, we adopted the Monte-Carlo 
method to generate a large number of simulated Lya 
forest spectra of type 1 bgQSOs that may be affected by 
(fg)QSOs, according to the statistical distributions of the 
HI column density, the Doppler parameter of the Lya ab- 
sorption lines, and the number evolution of the absorbers 
associated with these absorption lines. We extracted the 
distribution of the flux decrement (DA) or the optical 
depth near the fgQSOs from these spectra. We illus- 
trated that the effects of the fgQSOs on the bgQSO spec- 
tra depend on the fgQSO lifetime, the anisotropy in the 
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Fig. 10. — The expected transverse proximity effect of type 2 
QSOs at redshift ~ 4. Legends for panels (a) and (b) are the 
same a s thos e for Figure |9ji and d, respectively. See details in 
Section 14331 

fgQSO UV radiation, and the density enhancement in 
the fgQSO proximity region, given the strength of the 
UVB. According to these calculations, the PE of QSOs 
may be divided into the ionization proximity effect (IPE) 
and the density proximity effect (DPE). The absorption 
near a QSO may be enhanced due to the DPE but de- 
crease due to the IPE. Depending on detailed properties 
of QSOs, the combination of DPE and IPE can result in 
either an excess or a decrease of absorption in the near 
zone of a fgQSO. 

The LOSPE of a QSO with known luminosity is af- 
fected by the density enhancement in the QSO near zone 
but irrelevant to the QSO lifetime and the anisotropy 
in the QSO UV radiation, which can thus provide con- 
straint on the effective density enhancement in the near 
zone of the QSO. Based on the measurements of LOSPE 
for a low- redshift low- luminosity QSO sample (i.e., z '-^ 2 
and L^, - 5 X lO^ierg s'^ Hz"!; see KT08) or a high- 
redshift high luminosity QSO s ample (i.e., z ^ 4 an d 
L^o 8 X lO^^erg s"^ Hz"!; see lGuimaraes et al.|[2007l ). 
we obtained the effective density enhancement in the 
near zones of the QSOs in these two samples by matching 
the numerical results to the observations, respectively. 
Assuming that the density enhancement in the near zones 
of the fgQSOs is the same as that obtained from the 
LOSPE, the TPE has been simulated for both type 1 
and type 2 fgQSOs with a wide range distribution of the 
fgQSO properties, i.e., the lifetime and the half opening 
angle of the associated tori. Our numerical simulations 
show that the response of the TPE to the change in the 
fgQSO lifetime (or the luminosity evolution or fluctua- 
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Fig. 11. — The expected difference in the transverse proximity 
effect of type 1 fgQSOs at redshift ~ 4 from that of type 2 fgQSOs. 
The optical depth distribution near type 1 fgQSOs is denoted by 
r(I), and that near type 2 f gQS Os denoted by t(II). Legend s are 
the same as those for Figure [TOl See details in Section 14.2.31 

tion on timescale of 10^ — 10^ yr) and that to the change 
of the half opening angle of torus are different. The dif- 
ferences in the resulted TPE for different settings of the 
QSO lifetime are large for the high-luminosity samples, 
but they are relatively small for the low-luminosity sam- 
ple. These results suggest that the density enhancement 
near the fgQSOs, the fgQSO lifetime and the half opening 
angle of the tori associated with the fgQSOs can be simul- 
taneously constrained by the LOSPE and the TPE com- 
bined together of a sample with several hundreds bright 
QSOs. If a significant number of type 2 QSOs can be re- 
vealed by future observations, our simulations show that 
the TPE of type 2 QSOs can be significantly different 
from that of type 1 QSOs, and an observational search 
for the contrast will further help to distinguish the dif- 
ferent effects due to QSO properties or environment and 
improve the constraints. 

The lifetime and the opening angle of the torus are 
of fundamental importance for our understanding of the 
QSO physics and the growth of massive black holes 
(MBHs). The lifetime of QSOs is an important pa- 
rameter characterizing the luminosity evolution of QSOs 
and thus the assembly history of MBHs (including the 
mass and the spin evolution). A number of indepen- 
dent arguments suggest that the lifetime of bright QSOs 
is on the order of lO'' — lO^y r fe.g.. iHaiman fc Hull 
[Ml"; 'Martini & Weinberg '200T: Yu & Tremaine" '20021; 
^tcidcl et al. 2002; Yu & Lu 2004; Marconi ct al. 200l 
IYu fc Lul l2008: Shen 2009; Shan kar et all 120101: also see 
a review bv , Martini 2004). However, in principle, MBH 
growth may occur mainly through a long period contin- 
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uous accretion or alternatively many short-period accre- 
tion episodes; and those numbers are mainly obtained 
from demography of QSOs and they may only represent 
the net lifetime that the QSO luminosity in some wave- 
length range larger than a threshold, n ot the time pe - 
riod for each accretion episode (e.g., see iMartinil l2004l ) . 
As we have demonstrated that the TPE offers a way, 
if not the only w ay, to constrain t he episodic lifetime 
of QSOs (see also IAdelbergeH[2004D . which should pro- 
vide considerable insight into our understanding of the 
growth history of MBHs. The presence of a torus sur- 
rounding each QSO is proposed to be the underlying 
reason that leads to the classification of type 1 QSOs 
and type 2 QSOs. The half opening angle of the torus 
is an important parameter that determines the ratio of 
type 2 to type 1 QSOs. As demonstrated in this paper, 
the TPE, together with the LOSPE, can provide con- 
straints on this angle, which is independent of those pre- 
vious ways mainly thr ough observational ratios of type 
2 to t ype 1 QSOs (e . g..lTreister fc Urr^ [2nM iGiUi et al.l 



I to t ype 1 MSUs (e . g..rireister & Urrvli:^uua luun et al.l 
I2007t iTreister et al.ll2010[ ). As an alternative potential 



tool to constrain the opening angle of the torus, the QSO 
PE may be helpful in our understanding of the unifica- 
tion picture of different types of the AGN population and 
its underlying physics. 

KT08 measured both the LOSPE and the TPE for 
a low-redshift low-luminosity QSO sample (z ~ 2 and 
i^ufi ~ 5 X 10^°erg s~^ Hz~^), and their measure- 
ment on the TPE suggested that the QSO lifetime is 
< 10^ yr. That is apparently in contradiction with 
the detections of the TPE through other observations 
on He II or metal lines, which require a longer QSO 
lifetime tq > (2 — 3) x lO'' yr Ce.g.. iJakobsen et al.l 



200l IWorseck fc WisotzkH [20061: IWorseck et all [20071: 
Goncalves et al.ll2008[ ). T"hrough our Monte-Carlo simu- 



lations, we have shown that the apparent contradiction 
can be solved after considering a combination of the ef- 
fects due to density enhancement in the near zone of the 
fgQSOs and the obscuration of the tori associated with 

the fgQSOs^ 

Guiniar aes et al.l ()2007[ ) measured the LOSPE for a 
high-redshift high-luminosity sample (z ~ 4 and L^^ 



8 X lO'^^erg s^^ Hz^^), and their LOSPE measurement 
suggests that the density in the near zones of those 
QSOs is enhanced. In order to reproduce their mea- 
surement on the LOSPE, a density enhancement is also 
requi red in our Mont e-Carlo simulations. Compared 
with iGuimaraes et al.l (2007 ) estimates on the density 
enhancement, the required density enhancement here is 
smaller by a factor of 1.5 — 5 at a distance of < 
10 Mpc-1 Mpc but similar at a distance of > 15 Mpc. 
The main reason for this difference is that the Pois- 
son variance of absorption lines in the near zones of 
QSOs, which is included in our Monte-Carlo simulations, 
leads to an averaged absorption that is larger than a 
simple estimation obt ained without consider ing of the 
Poisson variance (see IGuimaraes et al.l l2007l ). There- 
fore, a less significant density enhanc ement is required 
here c ompared with that obtained by IGuimaraes et all 
(l2007f) (see also discu ssions on the Poisson variance in 
iDairAglio eraII[20Ci8h . 

We conclude that the current measurements on the 
LOSPE and the TPE are consistent with that (1) 
the density is significantly enhanced in the vicinity 
of fgQSOs, (2) the fgQSO lifetime is longer than 
a few 10^ yr and a much shorter lifetime (i.e., 
< 10^ yr) is excluded or at least is not required, 
and (3) the half opening angle of the tori associ- 
ated with fgQSOs is ~ 60°, which are consistent 
with other independent estimate s (e.g.. IJakobsen et al.l 
[20031 Worscck & Wiso t^ [200l IWorseck et al.l 120071: 
'Concalves et al. 2008; Treister fc Urrvll2006t iGilh et al.l 
^007; Treister et al. 201^^ Future observations on the 
PE of several hundreds of bright type 1 and type 2 
(fg)QSOs will potentially set accurate constraints simul- 
taneously on these QSO properties. 
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APPENDIX 

INTERSECTIONS OF THE LIGHT CONE CONFINED BY THE TORUS AND THE SURFACE CONFINED BY THE 

TIME-DELAY EFFECT 

The light cone confined by the dusty torus associated with the central engine of a QSO is given by 



+ (Y cos *o - R\\ sin *o) 
-(r sin^fo + i?|| cos*o)^tan^ Oq = 0, 



(Al) 



where (X, Y, i?||) are the three-dimensional Cartesian coordinates with the QSO being the origin, Qq is the half opening 
angle of the torus, \l/o is the offset of the torus axis from the LOS of the fgQSO (see Figure [1]). The cylinder with a 
radius of R\ is defined as 



(A2) 



where X — R± sin$, Y — R± cos*!", and $ e [0°, 360°) is the rotation angle around the LOS. Figure [T] shows the cases 
for $ = 0°. For any given $, the intersections of the light cone with the cylinder are the solutions of equations I All and 
\M[ If $ = 0°, the two solutions are e = -i?_L/tan(eo - *o) and i?||_F = i?x/tan(eo -I- *o) for type 1 QSOs (the 
left panel in Figure[T]) or e = i?_L/tan(8o -t- ^'o) and i?||^F = — i?_L/tan(6o — ^o) for type 2 QSOs (the right panel 
in Figure [Ij, respectively. For the case of $ = 90°, the two solutions are R\\,e = ^R±/{\ cos^t'olV tan^ Qq — tan^ 4*0) 
and i?||,F — ^R\\.E if |tanOo| > |tan^oh and there is no solution to the above two equations, otherwise. For those 
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cases without solutions to the above two equations, there is no the PE on the bgQSO spectrum due to the dusty torus 
associated with the type 2 fgQSO. 
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